Curvaton reheating allows TeV Hubble scale in NO inflation 
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Curvaton reheating is studied in non-oscillatory (NO) models of inflation, with the aim to ob- 
tain bounds on the parameters of curvaton models and find out whether low scale inflation can be 
attained. Using a minimal curvaton model, it is found that the allowed parameter space is consid- 
erably larger than in the case of the usual oscillatory inflation models. In particular, inflation with 
Hubble scale as low as I TeV is comfortably allowed. 



Inflation is to date the most compelling solution for 
the horizon and flatness problems of big bang cosmology. 
Recently, observational data at high precision have con- 
firmed that structure formation and the CMB anisotropy 
are due to the existence of superhorizon curvature pertur- 
bations, which need an acausal mechanism to be created 
and, therefore, arc further evidence of inflation. 

The majority of inflationary models suggest that infla- 
tion took place at energy comparable to that of grand 
unification. This is because, in most inflationary mod- 
els, the curvature perturbations are due to the quantum 
fluctuations of the inflaton field; the field controlling the 
dynamics of inflation. In this case the inflationary energy 
scale is [l[ 



V» 1/4 = 0.027ei /4 m P 



(1) 



where mp = 2.4 x 10 18 GeV is the reduced Planck mass 
and e = —H/H 2 is one of the so-called slow- roll param- 
eters, where H = a/a with a being the scale factor and 
the dot denoting derivative with respect to time. In the 
above, '*' denotes the epoch when the cosmological scales 
exit the horizon during inflation. For inflation to occur 
one needs e < 1, which is equivalent to a > 0. In slow- 
roll models of inflation e< 1, which means H ~ constant. 
However, slow-roll inflation is not exactly de Sitter ex- 
pansion, because it has to end for the hot big bang to be- 
gin. Therefore, albeit small, in most models e is not tiny. 
In fact, it turns out that, typically, e* ~ 1/iV*, where 
./V* is the remaining number of inflation e-folds when the 
cosmological scales exit the horizon. iV* ranges between 
40-65 for inflation scales between grand unification and 
Big Bang Nucleosynthesis (BBN). Hence, from Eq. (fl]) 

we see that v}^ 4 ~ 10 15 ~ 16 GeV. Exceptions to this rule 
exist, where e is exponentially suppressed due to appro- 
priate features of the scalar potential (e.g. see Ref. 
but they are typically fine-tuned (e.g. see Ref. Q). 

However, in recent years, advances in string theory 
have turned the attention at progressively lower energy 
scales. In particular, large extra dimensions suggest that 
the fundamental scale may by much smaller than mp 



and the scale of grand unification, possibly near the elec- 
troweak scale instead (e.g. see 3). This means that the 
hot big bang cannot be extended up to grand unifica- 
tion energies and so inflation must take place at lower 
scale. Moreover, in the search for the inflaton field, the 
most promising candidates are either flat directions in 
supcrsymmetric theories or string moduli fields, which 
correspond to the size and shape of the compactified 
extra dimensions. However, in the MSSM and its ex- 
tensions the corresponding scales are typically much less 
than the grand unified scale. Similarly, for string mod- 
uli one expects the scalar potential to vary significantly 
(i.e. 6V/V ~ 1) over distances ~ mp in field space, 
which means that the typical scale for modular inflation 



is K 1/4 ~ v / m 3/2 m P - 10 10 - 5 GeV (i.e. i?* - 1 TeV) for 
moduli with masses of order m^/2 ~ 1 TeV, such as string 
axions. Recently, many possibilities for inflation in string 
theory have arisen in the context of the string landscape 
For example, metastable supersymmetric vacua @ 
may account for the vacuum density of inflation. This 
sets the value of U* at the supersymmetry breaking scale 
which corresponds to the above mentioned scale or lower. 
This is why there has been growing interest in achieving 
inflation at low energy scales. 

One way to liberate inflation from the constraint in 
Eq. |T]) is to consider that the curvature perturbations 
generated by inflation are due to the quantum fluctu- 
ations of a field other than the inflaton, in which case 
Eq. |T]) turns into an upper bound Q. This, so called cur- 
vaton field, docs not influence the dynamics of inflation 
but becomes important after inflation has ended, when 
it imprints its curvature perturbation onto the Universe 
§. Under this hypothesis, it is possible to relax this 
constraint substantially when the curvaton model is spe- 
cially designed to do so [1, [T3|. However, generically, for 
a minimal curvaton model, there is still a lower bound 

on the inflationary scale of order v}^ 4 > 10 12 GeV (i.e. 

> 10 6 GeV) pT| . when considering inflation, which 
ends through oscillations of the inflaton field, as usually 
assumed. 

In this letter we show that, for a minimal curvaton 
model, the inflation scale can be substantially lower when 
considering a non-oscillatory (NO) model of inflation. In 
such a model the inflaton corresponds to a flat direc- 
tion with runaway behaviour, typical, for example, of the 
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scalar potential of the so-called geometric string moduli 
fields. In such models after inflation ends the Universe 
remains dominated by the inflaton field which is now fast- 
rolling down its runaway scalar potential [l2l ]. Since the 
inflaton does not oscillate and does not decay after the 
end of inflation, reheating occurs by other means. 1 Re- 
cently, a curvaton field has been assumed to be respon- 
sible for reheating in NO inflation models [HI, [lil, Il5j|. 
However, up to now the parameter space for the infla- 
tionary scale (in particular low-scale inflation) has not 
been studied. 

We begin by outlining the particulars of NO infla- 
tion. In NO models, after the end of inflation, the in- 
flaton field becomes dominated by its kinetic density 
P<i> ~ Pkin = \<\? [12]. The equation of motion for the 
field becomes cf> + 3H(f> ~ 0, which is oblivious to the 
scalar potential and suggests that oc a~ 6 . This gives 
rise to a period of so-called kination, where a oc 
Because the dilution of the inflaton's density is so dras- 
tic, the Universe can be reheated by the thermal bath 
due to the decay products of a curvaton field, which also 
accounts for the density perturbations in the Universe 
[lH . We should stress here that our treatment below is 
independent of the NO inflation model since the scalar 
potential for the inflaton sector is negligible during the 
evolution of the curvaton field. 

We consider a minimal curvaton scenario where the 
curvaton field a is characterised by the scalar potential 



V{a) 



f 



2 2 

to a . 



(2) 



In order to act as a curvaton a must undergo parti- 
cle production during inflation. Hence, it needs to be 
effectively masslcss. This implies the bound m < 
An even tighter bound, however, is due to spectral in- 
dex considerations. Indeed , for the curvaton we have 
n s - 1 ~ It] = | (to./H*) 2 , where j] = rrip^^ [|| with 
V* ~ 3H 2 mp being the density during inflation. Since 
observations do not favour a blue spectrum, we require 
that n s is at most n s ~ 1 .00, which is still marginally 
acceptable. 2 This means that to must satisfy the bound 



to < 0.1 if. 



(3) 



We assume that the inflaton's contribution to the cur- 
vature perturbation is negligible, as is the norm under 
the curvaton hypothesis [8j. In any case, it is difficult 
to obtain a sizable curvature perturbation in low scale 
inflation (See, however, Ref. @). Since the Universe 
is reheated by the decay products of the curvaton which 



1 NO models have been particularly useful for quintessential infla- 
tion (e.g. see Ref. Il3t 1 14l and references therein) because the 
inflaton field can survive until today and be responsible for the 
dark energy at present. 

2 Smaller values of n s can be attained if e is not tiny, in which case 
n s - 1 ~ 2(r? - e) [|. 



dominate the scalar field background, the observed cur- 
vature perturbation £ = 4.8 x 10 -5 is entirely due to the 
curvaton, i.e. 
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where we used that, before its decay (denoted as 'dec'), 
the curvaton undergoes quasi-harmonic oscillations in the 
scalar potential of Eq. ([2]), so that p a = 2V(a) oc a~ 3 
[ijj, where V(a) is the average value of the potential den- 
sity over many oscillations. Since the oscillating curvaton 
satisfies the same linear equation of motion as its oscillat- 
ing perturbation, (namely a + 3H& + m 2 a = 0) the frac- 
tional perturbation 5a /a remains constant and equal to 
its value at the onset of the oscillations, denoted as 'osc'. 
Before the onset of the oscillations, both the value of the 
field and its perturbation are frozen so the fractional per- 
turbation remains constant again and equal to its value 
at horizon exit, denoted by '*'. As is the case with parti- 
cle production of scalar fields, the value of the curvaton 
perturbation at horizon exit is determined by the Hawk- 
ing temperature: <5er» = H*/2ir. Note that Eq. Q im- 
plies that a* 3> i?* , which guarantees that the curvature 
perturbations are predominantly Gaussian Q. 

After the end of inflation, the inflaton becomes ki- 
netically dominated, while the curvaton remains frozen 
because to <C -ff* and its motion is overdamped. This 
means that the density parameter of the curvaton 
SI = Pa/ P scales as 17 oc a 6 until H(t) drops down to 
H ~ to, when the curvaton unfreezes and begins its 
quasi-harmonic oscillations. At this time we have 



( a * 



(5) 



After the onset of the oscillations, the curvaton density 
decreases as Pa oc a~ 3 so that f2 oc a 3 . 

Now, there are two possibilities, depending on whether 
the curvaton decays before its density dominates the Uni- 
verse or afterwards. Let us assume firstly that the cur- 
vaton decays before domination. In this case, using that 
during kination a oc H^ 1 ^ 3 , one readily obtains that the 
curvaton density parameter at the time of curvaton decay 
is 

2 



TO / 67* 

r \mp 



(6) 



where T is the curvaton decay rate. The decay products 
of the curvaton redshift as relativistic matter p 7 oc a~ 4 , 
which means that their density parameter grows as 
f2 7 oc a 2 . Reheating is achieved when this density param- 
eter reaches unity and the Universe becomes dominated 
by the thermal bath of the curvaton decay products. Us- 
ing that i/dom ~ 2~reh/ m -P' ^ 1S eas y to find 



3 



where we also employed Eq. §4§ anc l we denoted with 
'dom' the epoch when the inflaton's density becomes sub- 
dominant. From the above and using also that T > -f/dom 
and T ro h > Tbbn we obtain the bound 



H 



>(10C 2 ) 1 / 3 (T| BN m F ) 1 / 3 ~10GeV 



(8) 



where we also considered Eq. ([3]) and demanded that re- 
heating occurs before BBN; the latter corresponding to 
temperature Tbbn ~ 1 MeV. Thus, we see that inflation 
with H* ~ 1 TeV is possible to accommodate with a min- 
imal curvaton model. 

Let us consider, now, the case when the curvaton de- 
cays after it dominates the Universe. In this case the 
oscillating curvaton dominates when 



mp 



(9) 



Reheating is, now, achieved at the decay of the curvaton, 
so that 



Trch ~ V mpY . 



(10) 



Using that T < i/dom and also T re h > Tbbn, we arrive 
once more at the bound in Eq. ([5]), where we also con- 
sidered Eq. 

Additional bounds on the parameters are obtained as 
follows. For the decay rate of the curvaton we can write 
r = g 2 m, where g is the coupling of the curvaton to its 
decay products. The expected range for this coupling is: 



m <r <r 1 

<; <7 <; 1 

mp 



(11) 



where the lower bound is due to gravitational decay. 

Assume at first that the curvaton decays before domi- 
nation. Then, from Eq. (J7J), we have T re h oc T -1 / 4 . Con- 
sequently, combining the lower bound in Eq. with 
the requirement that T ro h > Tbbn and Eq. ([7]), we ob- 
tain the bound 



H*> ((Ti^mp) 1 / 3 ~ 1 GeV, 

which is weaker than the bound in Eq. 
that T rc h > Tbbn, Eq. ([7]) results in 



.9 < 



mpm 

1 BBN 



/ / 

(rn P 



< io 41 c 3 



mp 



(12) 
Also, using 



(13) 



where, in the last inequality, we used Eq. Com- 
paring this bound with the lower bound in Eq. ([lip it 
can be easily verified (with the help of Eq. (fT2"| ) that 
there is always parameter space for g. Using also that 
g < 1 < 0.1 H*/m according to Eqs. and (fTTj) , it can 
be shown that the bound in Eq. ([T3|) is tighter that 
the upper bound in Eq. (JTTJ) only for H* < 10 TeV. Fi- 
nally, since in the case of decay before domination we 
have T > f/dom ~ T 2 eh /m P , employing Eq. the upper 
bound in Eq. (fTTj) suggests 



< 



mp 



H*<(m P ~ 10 i4 GeV, 



(14) 



where we also used Eq. (|4|). 

Now, let us consider the case when the curvaton de- 
cays after domination, where T re h is given by Eq. (|10[) . 
Combining the requirement that T rc h > Tbbn with the 
upper bound in Eq. (|lip and using also Eq. one 
finds the trivial bound: ff» > 10 Tg BN /mp, which is way 
weaker than the bound in Eq. ([5]). The corresponding 
lower bound on g is 



1/2 



Tbbn 
mp 



(15) 



which is tighter than the lower bound in Eq. (jlljl if 
m < (^BBN TO f ) 1 ^ 3 - Finally, combining the condition 
r < i?dom with the lower bound in Eq. (jllj) one ends 
up again with the bound in Eq. p4[) . 

In summary, the allowed range for the inflationary 
Hubble scale is 

10 GeV < H* < 10 14 GeV, (16) 
while, for the decay coupling, the allowed range is 



Tbbn m 
\ y/mpm ' mp 



mpm 



( — Y 



> ~ 9 ^min< 1, ■ . 

t bbn \mpj 



(.17) 



Another set of bounds is due to the possible overpro- 
duction of gravitational waves (GWs) due to inflation. If 
inflation is followed by a phase whose equation of state 
is stiffer than radiation, the spectrum of relic gravitons 
features a spike (the slope grows with the frequency) for 
modes re-entering the horizon during the stiff phase [l8| . 
Since in NO models there is a phase of kination right 
after inflation, high frequency gravitons re-entering the 
horizon during kination may disrupt BBN by increasing 
H . To avoid this, we require [H, E3 



r k, 

I = h 2 Q GW (k)dliik < 2 x 10" 

■' k-RRXI 



(18) 



where ilow(fc) is the GW density fraction with physical 
momentum k and h = 0.73 is the Hubble constant Ho in 
units of 100 km/sec/Mpc. Using the spectrum Ogw(^) 
as computed in [181 ]. the above constraint can be written 
as fil 



1 Hi 



1 1 



H 



dom 



2/3 



(19) 



where acw — 0.1 is the GW generation efficiency dur- 
ing inflation and f2 7 (fco) = 2.6 x 10~ 5 h~ 2 is the density 
fraction of radiation at present on horizon scales. 

If the curvaton field decays before domination, using 
Hdom ~ T t 2 oh /mp and Eq. JjJ, the bound in Eq. (JTSj) be- 
comes 



^(-L)" 3 <24C 
m \H* 



10 



10 



(20) 
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If the curvaton decays after domination, using Eqs. ([9]) 
and pop , the bound in Eq. (Tg) becomes 



The above bounds may truncate further the ranges in 
Eqs. (32) and (31). 

Let us quantify the above with a couple of spe- 
cific examples. Firstly, let us choose H * ~ 1 TeV and 
m ~ 100 GeV so that the bound in Eq. (3) is saturated. 
In this case Eq. (4) suggests that er* ~ 10 s GeV. Using 
this, Eq. (3D becomes 10" 13 < g < 10~ 4 . Using Eq. (5) 
we find i/dom/r ~ 10 _20 /<7 2 , which means that the cur- 
vaton decays before domination if 10~ 10 < g < 10~ 4 , 
whereas it decays after domination if 10~ 13 < g < 10 -10 . 
In the former case Eq. (7) gives T rch ~ ,g~ 1/2 10~ 5 GeV, 
while in the latter case Eq. (flO) gives T roh ~ glO 10 GeV. 
Hence, the allowed range for the reheating temperature 
is: 1 MeV < T lch < 1 GeV. Since T = g 2 m, it is straight- 
forward to check that the bounds in Eqs. (p?0") and (pL) 
are satisfied in the range shown in Eq. (17)) . 

Now, let us choose cr* ~ mp and m ~ 100 GeV, cor- 
responding to a string axion as curvaton. In this 
case Eq. (?) suggests that ~ 10 14 GeV, which sat- 
urates the bound in Eq. (14) and corresponds to 
inflation at the grand unified scale. Using this, 
Eq. (37) becomes 10~ 13 <g<l. Then, Eq. © gives 
fldom/r ~ .9~ 2 ^ 1- This means that the curvaton has 
to decay at or after domination. Thus, from Eq. (TO) 
we have T ro h ~ glO 10 GeV, which results in the range: 
1 MeV < T re h < 10 10 GeV. This time, however, the GW 
constraint turns out to be much tighter. Because the 
curvaton decays at or after domination, we need to use 
Eq. pT) . With the chosen values, the bound can only be 
satisfied with g ~ 1, which saturates the upper bound 
in Eq. (TT). Thus, the reheating temperature has to 
be T re h ~ 10 10 GeV, which seriously challenges gravitino 
constraints. 

In conclusion, we have studied in detail the parameter 
space for the inflation scale in NO models with a min- 
imal curvaton scenario being responsible for reheating 
the Universe as well as for the curvature perturbations. 
We have shown that low-scale inflation with as low 
as 10 GeV is feasible in this case, regardless of whether 
the curvaton decays before or after dominating the Uni- 
verse. Our results are independent of the particular form 
of the NO inflation model because, during the curvaton 
evolution, the inflaton sector is oblivious of the scalar po- 
tential, since it is dominated by the kinetic density of the 
inflaton. 
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